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Abstract 

We propose a Wigner function based parameter that can be used as an indicator of quantum 
chaos. This parameter is defined as "entropy" from the time-dependence of " non-classicallity" 



proposed in [lj. We perform our considerations for the system of damped nonlinear (Kerr-like) 
oscillator excited by a series of ultra-short external pulses. 
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I. INTRODUCTION 



The systems demonstrating chaotic behaviour in their dynamics are still the subject 
of much attention. Especially, there has been a great deal of interest in exploring classical 
dynamics of chaotic systems. For such systems certain methods allowing distinction between 
the regions of regular dynamics from those of chaotic one are well developed and widely 
used. It is still of special importance to find the adequate methods for analysis of the 
systems whose quantum dynamics exhibits chaotic behaviour. It is of interest not only from 
the cognitive point of view but is also for the development of quantum information theory 
methods. It is known that chaotic behaviour in the dynamics of a quantum system should 
destroy the entanglement (which is the essential point in the quantum computing) between 
the quantum states, but surprisingly in some cases quantum chaos can even enhance the 
stability of quantum computation 2|. 

Many attempts to address the problem of quantum chaos have already been made. For 
example, it is known that there is a correspondence between the statistics of eigenvalues 
and eigenvectors of quantized classically chaotic systems and the canonical ensembles of the 
random matrix theory {3], [5], [f], Q, g[ Is]. The distances between successive eigenvalues 
of quantized chaotic system have the same probability distribution as the successive dis- 
tances between the eigenvalues of random matrix. There is also a method based on the 



fidelity decay, indicating the existence of chaos in dynamics of quantum systems [10l. I 111. Il2|. 
The time-evolution of the fidelity between the state evolving under the unitary mapping 
procedure and the same initial state evolving under the same map but subjected to some 
additional tiny perturbations, in the chaotic region exhibits an exponential decay [12j]. It 
means that there is a significant difference between these two quantum states and that the 
dynamics of the quantum system is sensitive to the initial conditions. This sensitivity is a 
characteristic feature of the chaotic dynamics. 

In the present article we would like to apply the quantum parameter based on the Wigner 
function (and connected to its negative volume), introduced in [lj] and referred to as the "non- 
classicallity parameter" , to the problem of a quantum counterpart of the classically chaotic 
system dynamics. We introduce an entropic measure based on this parameter and show it to 
be sensitive to changes in the regions where the classical counterpart of the system described 
is chaotic. These changes indicate that the final state of the quantum system in some regions 
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is also very sensitive to the initial condition variations and hence, very tiny perturbation of 
the initial conditions lead to considerable changes in the final quantum state. To investigate 
such phenomena we would apply the recurrence plots analysis. It will allow us to confirm 
our statements about the chaotic behaviour (in the classical sense) of the strictly quantum 
parameter defined from the Wigner function of the system described. 



II. THE MODEL 



We consider a system composed of an anharmonic oscillator driven by a series of ultra- 
short external pulses. Various aspects of such a system's dynamics have been discussed in 
many papers (for example see Q, Q] and the referrences quoted therein) . It is known 
that depending on the values of the parameters used, the system can exhibit regular or 
chaotic dynamics The problem of the classical dynamics and that of the quantum 

counterpart of this system is also discussed in that paper. 

The system is supposed to be initially in the vacuum state. In the interaction picture, the 
time evolution of the system under consideration is governed by the following Hamiltonian 
(we use units ft = 1): 

H = H NL + H K , (1) 

where Hnl describes the evolution between subsequent pulses and Hk - the evolution 
initiated by an ultra-short pulse: 



H NL = ±(a + ) 2 a\ (2) 



2 (° ) a 

H K = e (a+ + a) 5 ( f ~ kT ) ■ ( 3 ) 

h=l 

Operators a + and a are those of photon creation and annihilation, respectively; \ describes 
the non-linearity of the oscillator (for Kerr medium it is the third order susceptibility) and 
in further considerations it is set to 1; e is the strength of the external pulse - nonlinear 
oscillator interaction, and T is the time between two subsequent pulses. Assuming that 
time T exceeds significantly the the duration of single pulse, we can model the series of 
these ultra-short pulses by a series of Dirac-delta functions. 

In the present considerations we deal with the damping case only. It seems to be a more 
realistic situation to let an oscillator to interact with the enivronment than simply deal with 
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a no-damping case. In real physical situations, the influence of the environment (represented 
by the damping process) may significantly change the system's dynamics. For the system 
considered the influence of the damping process makes the system's dynamics regular for 
a wider range of external excitation strength than in the no-damping case. As the system 
is damped in the chaotic region fewer states are involved in the system's dynamics, which 
significantly simplifies numerical calculations. When we do not include the damping process, 
in the chaotic region many more states have to be considered and additionally, the chaotic 
behaviour in the quantum system considered is visible for longer times. If we look closer 
at the numerical complexity of the problem it becomes clear that for the cases without 
damping it would be more suitable to analyse for example a fidelity decay as a signture of 
chaotic dynamics in a quantum system rather than to deal with a Wigner function based 
parameter. When dealing with a Wigner-based parameter (which has to be calculated and 
analysed after each external kick) it would be more desirable to work with smaller basis and 
to have an opportunity to observe chaotic behaviour of the system in shorter times as well. 

For the system considered the density matrix approach is necessary and we use the 
formalism proposed in [uj and used for example in 16]. The time evolution of the density 
matrix is governed by the following master equation: 

■J^ = — [(a + ) 2 a 2 , p] + ^ (2apd + — a + ap — pa + a) , (4) 

that can be solved analytically (for example see [l?}]) or numerically. In fact, equation 
flU) describes the free evolution of the damped system between two subsequent pulses. To 
include the interaction with the external field, one should apply to the density matrix (after 
its free evolution during the time T) the operator Uk of the form: 

U K = e - i <* + +") . (5) 

Consequently, the whole dynamics will be described successively by the evolution (during 
the time T) according to the master equation (j3J) and the "kicked" operator (jSJ). The initial 
state of the system considered is the ground state pit = 0) = |0)(0|. 

It is known that the classical kicked Kerr oscillator, depending on the values of the pa- 
rameters used, can exhibit regions with classically regular or chaotic behaviour. Appropriate 
bifurcation diagrams of the mean classical energy as a function of the external coupling pa- 



rameter e, for weak and strong damping cases have been presented in 



161 ] . The procedure 



used for obtaining such diagrams first requires the explicit form of the equation of motion for 
the annihilation operator for the time between the subsequent external pulses. The solution 
of this equation can be expressed in a simple analytical form, because at first we neglect 
the damping process, and consequently, the number of photons h = a + a is conserved during 
the evolution. Then, we include the influence of the external pulse (described by the action 
of the operator (jSJ) which in fact is a shift operator. The final point of the procedure is to 
replace all the operators (a + , a) by the complex numbers (a*, a) and the damping rate 7 
can be introduced at this point. Finally, the equation for a has the form: 

a k+1 = (a k -te)e-<^-^-^ T , (6) 

and the formula for the classical mean energy | a \ 2 needed for the construction of the bifurca- 
tion diagram can be obtained easily. A classical state which would be a classical counterpart 
of the initial system state (the vacuum state) can be simulated as in [l^]. To introduce the 
initial condition we have created an ensemble of classical trajectories whose starting points 
had been randomly chosen from a circle of radius 0.5 and centred at a = and after that 
the average trajectory has been treated as a classical counterpart of a quantum one. 

In the present paper we shall concentrate on the case when damping constant 7 = 0.1. 
The bifurcation diagram for a classical mean energy (Fig.l) shows that for a wide range of 
the external kicks strength e <~ 0.9, one can observe the regular dynamics of the classical 
oscillator. The chaotic behaviour appears when e exceeds the value of ~ 0.9. Additionally, 
the bifurcation (star ting from e ~ 0.7) appears in the diagram. This situation corresponds 
to that discussed in lfjj. Moreover, one should remember that the regular and chaotic be- 
haviour visible in the bifurcation diagram concerns the classical dynamics of the anharmonic 
oscillator. 

In the present paper we will try to answer the question about the regularity or chaocity 
in the dynamics of the quantum system whose classical counterpart behaves regularly or 
chaotic. In other words we will try to make use of a quantum parameter (characterising the 
quantum system) to identify the regions of regular quantum dynamics and those of a quan- 
tum chaotic nature. We have chosen a quantum parameter (reffered to as "nonclassicality") 
based on the negativity of the Wigner function [l| for the analysis of quantum dynamics 
of the kicked nonlinear oscillator. As claimed in this parameter is connected to the 
quantum character of the states analysed and consequently, we shall check its usefulness as 
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an indicator of quantum chaotic behaviour. 



III. THE EVOLUTION OF THE NON-CLASSIC ALLITY INDICATOR 

This paper is devoted to the applicability of the "non-classicallity" parameter as an 
indicator of quantum chaos. This parameter has been introduced in |l| and is related to 
the negative volume of the Wigner function. It is known that the coherent states minimise 
the uncertainty principle and in this sense they are considered as classical states. When 
we talk about the "non-classicallity" of the states we understand this "non-classicallity" as 
the deviation from the coherent states. It is known that the Wigner function is a quasi- 
probability function that represents a quantum state in the coherent state basis. The Wigner 
function is joined with the symmetrically ordered bosonic operators. For a state described 



by the density matrix, p is defined as [181 ]: 



W{a) = 2Tr D _1 (a)pD(a)P , (7) 

where D(a) = exp[aa + — a* a] is the displacement operator and P = exp[iTra + a] is the 
parity operator. The Wigner function can take both positive and negative values. The 
non-classicallity indicator proposed and discussed in [l| is defined as follows: 

5(p) = J [\W(a)\-W{a))da . (8) 
The value of 5 depends on the volume of the negative part of the Wigner function and is 



equal to zero 
non-negative 



'or the coherent and squeezed vacuum states for which the Wigner function is 



rg| . In this sense, the 5(p) parameter describes how much the quantum state 
considered differs from the coherent one and in consequence, form the "classical" state. It 
was shown by the authors of the 5{p) definition, that the value of this parameter for the 
Fock states \n) grows monotonically with n and consequently, the higher the value of n the 
greater the deviation of the state \n) from the coherent state. Moreover, the value of 5(p) 
does not depend on the degree of squeezing. 

We have analysed the time evolution of the parameter 5(p) defined as in ([HD, concentrating 
on the parameters describing our system (T, x, 7), used for analysis of the bifurcation 
diagram. From the whole range of external field - nonlinear system coupling strengths e we 
have chosen the values corresponding to the classically regular and chaotic regions. 



Moreover, in this paper we use S n (p) = 5(p)/n instead of S(p), which means that the non- 
classicallity is divided by the mean number of photons after each external kick - so we use 
the parameter which can be more convenient in further considerations when the number of 
photons grows considerably (regions of chaos). In the regions of regular dynamics it causes no 
significant changes in the values of S(p). Figures 2 and 3 present the time-evolution of S n (p) 
for the excitation strengths corresponding to the classically regular e = 0.2; 0.3; 0.4 - Fig.2a 
and chaotic e = 1.24 — Fig. 2b dynamics. Thus, we have found that whenever the classical 
system exhibits regular dynamics, the changes in S n (p) are regular, and depending on the 7 
value they may be slower or faster damped. The character of time-evolution does not change 
with increasing value of e. The regular damped oscillations (whose frequency depends on 
the value of e) occur - this behaviour can be seen in Fig. 2a, where we have plotted S(p) 
versus time for various values of external kick strength e. One should keep in mind that for 
such values of e the classical counterpart of our quantum system behaves regularly. We have 
plotted the time dependence of the mean energy of the classical system in that regular region 
- Fig.3a (for recollection - an average of randomly choosen 10 5 - in our case - trajectories 
starting from a circle of radius 0.5 centered at Re(a) = Im(a) = 0). We can see that indeed, 
the mean energy of the system (found for a sufficiently long period of time) does not vary 
significantly with changes in the external pulses strength. For comparison, when we plot 
the mean number of photons of the quantum system — Fig. 3b — (which would be a kind of 
analog of the mean energy of the classical system) we can observe a similar behaviour of this 
quantity. Small changes in the kicking strength cause the same small changes in the final 
mean number of photons (n) in the system. For the case of strong external pulses (e > 0.9 - 
chaotic region) the situation changes considerably. The changes in the values of the system's 
final energy are more sensitive to the changes in the external pumping strength. While for 
the case of weak excitations a change in e of ps 0.1 causes a change in the mean number 
of photons of PS 0.09 for e > 0.9 the same changes in the pump strength cause 10 x larger 
increase in the mean number of photons generated in the process. Moreover, the energy of 
the quantum system initially increases rapidly (as in its classical counterpart) for the values 
of e corresponding to the "deep chaos" case in the classical system. We can clearly see the 
similarities of the essential features of the two quantities discussed here (mean energy and 
mean number of photons). The values of the energy for the classical and quantum system's 
are not exactly the same but the character of their changes is similar in both (classical and 
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quantum) models. 

More interesting is the situation when the pulses are sufficiently strong to lead to chaotic 
behaviour in the classical counterpart of the system described. We have chosen the kick 
strength e = 1.24. It is seen from the bifurcation diagram that this value is high enough 
to put the classical system into a region of deep chaos. For comparison, we can analyse 
the behaviour of the quantum counterpart of the kicked oscillator for the same values of 
the other parameters as in the previously discussed case. Fig. 2b presents the evolution 
of S n (p) characterising the quantum system for the parameters which caused the classical 
counterpart of the system analysed to behave chaoticly. It is seen that there are no regular 
oscillations any more. Irregularities apparent in the time-evolution of S n (p) indicate that the 
final state of the quantum system considered cannot be well defined and changes significantly 
from one external pulse to another. Many states with various photon numbers are involved 
in the process and probability none of them is dominant. We can see that the parameter 
S n (p) indicates almost instantaneous irregular changes. The range of values reached by S n (p) 
indicates that the Wigner functions have noticeable negative parts (which is a measure of 
the quantumnes of the state). The negative fractions of W also change from one pulse to 
another. 



IV. RECURRENCE PLOTS ANALYSIS 



We have identified irregular changes in the quantum parameter S n (p) for the cases cor- 
responding to the situation when the classical counterpart of the system exhibits chaotic 
behaviour. We suppose that these irregularities are connected with the fact that the quan- 
tum system dynamics is also chaotic. To confirm our supposition we will analyse the time 
dependence of the " non-classicallity" parameter via the recurrence plots (RP) method, which 
can be easily applied for a short time series. 

The idea of recurrence plots comes from [20] and was developed and widely used in 
nonlinear data analysis. This method can be used to solve problems related to medicine 



[see for example 



21 



221]). economy [23|, |24j, geophysics [25|, |26(, astronomy [27|, |28( and 



many other fields, whenever there is a need to analyse a time varying signal coming from 
a nonlinear system. The basic idea of RP is to graphically present the times at which 
the system analysed recurs, or more precisely, the times at which the system's trajectory 
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presented in an appropriate phase space returns (with some approximation) to the same 
area. RP allows analysis of nonlinear signals from the systems whose trajectories generally 
exist in a many-dimensional phase space via a 2-dimensional plot. 

To construct RP we need to determine the binary matrix R. Its elements are defined as 

Rij = Q(t t hr Xj\\) i,j = l---,N Xi 6 3? , (9) 

where is the Heviside function, e t h r is the threshold parameter and 1 1 ■ 1 1 denotes the norm. 
This norm allows determination of the distance between two points. The trajectory of the 
system (obatined as time series of some system's parameter — for the case discussed here 
it is S n (p)) is first of all reconstructed in an apropriate phase space (whose dimension is 
determined before) at time t. Consequently, Xi — Xj measures the distance between the 
points of the so reconstructed trajectory. If two points fall inside the same region (sphere) 
according to e^ r , we would label them by 1, otherwise they are labelled as 0. In consequence, 
one obtains a square matrix with zeroes and ones, or graphically with black and white points. 
The matrix can be analysed using the Recurrence Quantification Analysis (RQA) intro- 



duced by Zbilut and Webber Jr |2jJ, [30j, based on examination of diagonal structures of RP. 



This method has been extended by Marwan 311] .who has proposed new measures of complex 



ity based on vertical structures in RP. In general, long diagonal lines in RP are characteristic 
of periodical orbits, whereas homogeneously distributed black points indicate white noise. 
On the other hand, the system exhibiting chaotic dynamics will produce shorter diagonal 
structures in RP and isolated black points among them. Moreover, the vertical structures 
suggest the existence of laminar states (states which do not change or change very slowly 



with time 31]) 



In this paper we have applied the RQA method for analysing the time varying "non- 
classicallity" parameter 5 n (p). We have concentrated on the regions where the classical 
counterpart of the system analysed behaves chaotically and we suppose that the quantum 
system exhibits chaotic dynamics as well. For the damping constant 7 = 0.1 we have 
obtained irregular changes in S n (p) (see Fig. 2b). RP would help us determine whether the 
system characterised by the time varying parameter S n (p) exhibits chaotic dynamics or not. 



To generate the RPs we have used a CRP toolbox by Marwan 32]. 

The first task is to reconstruct the phase space trajectory of the system using a time 
series of a measured quantity (in our case the "non-classicallity" of the system). It can be 



done using a time-delay method, which involves the use of the apropriately chosen minimal 
sufficient dimension for reconstructing the original system's trajectory. At first, we have 
estimated the optimal time delay r (using mutual information function) as r = 1. It has 
to be chosen in such a way that the linear dependencies between two subsequent vectors in 
a reconstructed phase space are reduced. It can be realized by finding the minimum of the 
mutual information function which describes the joint probability of finding the time series 



33|. 



value at the i-th interval and after a time r in the j-th interval 

Moreover, we have found the embedding dimension as d = 10. We have used for this 
purpose the false nearest neighbours function that counts the number of points in the nearest 
neighbourhood of a specified point. A point marked as a neighbour of another one in lower 
dimension may not belong to the neigbourhood of this point in a higher dimension phase 
space - such a point is called a false neighbour and the false nearest neighbours function 
allows choosing such a phase space dimension at which the number of false neighbours 
vanishes. 

For these parameters we have generated the RP presented in Fig. 4a. We can see a pattern 
that is characteristic of chaotic dynamics — some diagonal lines are visible. Their length 
distribution is plotted in Fig.4b. The maximum length of the diagonal line is 26 points. 
Although the lines of shorter lengths are dominant, there is a significant amount of lines 
with lengths larger than 2 points. Formation of diagonal lines is characteristic of chaotic 
dynamics of the parameter analysed and their lengths are related to the value of the largest 
Lyapunov exponent of a chaotic system 20|, |34j . From RQA we have obtained the average 
diagonal line length of 4 points. It is worth stressing that most of the points that are present 
in the RP form diagonal structures - the percentage of these points is ~ 92%. Therefore, 
we can say that the process of non-classical states generation initiated by external pulses 
(appearing as an increase in S n (p)) during the time of the evolution, for the excitations 
sufficiently strong to cause chaos in the classical counterpart of the system considered, is 
rather of deterministic nature. Stochastic processes are characterised by isolated points or 
very short lines (mainly composed of 2 points) and a low value of DET - contrary to the 
case discussed here. 

Therefore, we can conclude that the parameter S n (p) can be treated as an indicator of chaos 
generated in quantum systems. 



10 



V. THE ENTROPIC MEASURE OF CHANGES IN NON-CLASSIC ALLITY PA- 
RAMETER. 



As we are interested in a strictly quantum parameter that would allow us to determine 
whether the system is in the regular or chaotic region of its dynamics, we would use the def- 
inition of the non-classicallity parameter S n (p) and propose a quantity (having the physical 
sense of entropy) that measures the changes in 5 n (p). We can define this entropy in a few 
steps: 

1. We define the Fourier transform of the 5 n (p) as: 



2. Then we calculate a "power spectrum" V = |JF(cj)| 2 and normalise V to get Vn{<^). 

3. Finally, we define the entropic measure of changes in the non-classicallity indicator 
through the power spectrum as: 



In consequence, we obtain one value of entropy for the specific coupling strength e. The 
entropy value equal to zero would indicate that the non-classicallity parameter 5 n (p) does 
not change in the whole range of time considered. The increase in the value of E and its 
irregular changes for various values of coupling strength would indicate significant changes 
in the time dependence of the S n (p) parameter. From the calculations of the above defined 
entropy, we have found that for the regions of regular system dynamics the changes in the 
value of the entropy with e are smooth in character (Fig. 5) indicating regular oscillations 
of S n (p). Changes in the value of entropy arise from the fact that the frequency of 5 n (p) 
oscillations change with e and for higher values of e more oscillations of S n (p) would appear 
prior to their vanishing (see Fig. 2a). Additionally, with higher values of e more states are 
involved in the process. But as long as the changes in E are smooth, the system's dynamics 
is regular. When the system is in the deep chaotic region we see significant and irregular 
changes in the value of entropy. These irregularities (appearing for excitation strengths 
corresponding to the classical deep chaos) are visible in Fig.5 (from e w 1). Therefore, we 




(10) 




(11) 
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can see that the entropy E (that is a strictly quantum parameter) can indicate irregular 
changes in the dynamics of the quantum system described. These irregularities appear for 
the values of e leading to the chaotic behaviour of the classical counterpart of the quantum 
system described. We have already analysed these regions via the recurrence plots of the time 
varying 5 n (p) parameter, and we have identified the regions corresponding to the classical 
deep chaos as being chaotic also for the quantum system. Therefore, we can conclude that 
the regions of irregular changes in the entropic parameter proposed are those of chaotic 
nature, indicating the presence of quantum chaos phenomena. 



VI. CONCLUSIONS 



The quantum systems, whose classical counterparts exhibit the chaotic behaviour are 
still widely explored. Their analysis is expected to answer the question whether the chaotic 
dynamics of a classical system would persist in the dynamics of its quantum counterpart. 
It is known that for the classical systems there are well known methods allowing distinction 
between the regular and chaotic regions, whereas for the quantum models such methods are 
still being developed. Some of them have been already proposed and discussed, and we have 
mentioned them in the introduction of this article. 

In this paper we have proposed an entropic parameter E (of strictly quantum nature) 
based on the non-classicallity proposed in connected with the Wigner function. We 
have shown that this parameter used for the analysis of the quantum system which in the 
classical limit can exhibit chaotic dynamics, undergoes rapid and significant changes with 
coupling strength e, in the regions corresponding to classical chaos. As an example we have 
analysed a quantum Kerr-like oscillator externally kicked by the ultra-short pulses. It is well 
known that nonlinear quantum systems can be used as the components of systems exhibiting 
quantum entanglement and as such can be considered as a tool in quantum informatics. It 
is therefore of special interest to explore simple methods allowing the identification of the 
regular or chaotic regions in dynamics of such nonlinear quantum systems. The behaviour 
of the entropic measure E for such a system indicates that the final state of the quantum 
system in the regions corresponding to the deep classical chaos cannot be well established 
and changes irregularly from one pulse to another. The entropic parameter, apart from the 
sensitivity to tiny perturbations, reveals irregular dynamics even for short times. We have 
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analysed this short-time dependence of the non-classicality parameter (which is the basis of 
the entropic measure) using the RP method - a tool used in time-series analysis. We have 
concluded that the dynamics of changes in the non-classicality of the state generated in the 
process discussed for the regions corresponding to classical chaos can be treated as being 
chaotic as well. Therefore, we believe that the parameter E proposed in this paper could be 
a useful tool for investigation of systems exhibiting quantum chaotic phenomena. 
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FIG. 2: Time dependence of the " non-classicallity" parameter S n (p) for various values of e, at (a) 
e = 0.2; 0.3; 0.4; at (b) e = 1.24. Other parameters: T = ir and 7 = 0.1 




FIG. 3: Time dependence of the mean energy of the classical system (a) and the mean number of 
photons (b) for various values of e. The parameters: T and 7 are the same as in Fig. 2 
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FIG. 4: Recurrence plot (a) for the signal presented in Fig.2b and diagonal line distribution at (b). 
Time delay is equal to 1, embedding dimension 10. The values of 7 and T are the same as in Fig. 2 
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FIG. 5: Entropy of changes for S n (p) versus excitation strength e. We assume T = ir and 7 = 0.1. 



